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We propose that localized defects in magnetic semiconductors act as deep impurities and can be 
described by the Anderson model. Within this model, hybridization of d-orbitals and p-orbitals gives 
rise to a non-RKKY indirect exchange mechanism, when the localized d-electrons are exchanged 
through both conduction and valence bands. For semiconductors with indirect band gap the non- 
RKKY part of exchange integral is antiferromagnetic, which suppresses ferromagnetism. In case 
of direct band gap, this exchange mechanism can, under certain conditions, lead to enhancement 
of ferromagnetism. The indirect exchange intergral is much stronger than RKKY, and can be 
sufficiently long range. Thus, a potentially new class of high-temperature magnetic semiconductors 
emerges, where doped carriers are not necessary to mediate ferromagnetism. Curie temperatures in 
such magnetic semiconductors are determined mostly by the interaction between localized impurities, 
not Zener mechanism. This effect could also be responsible for unusually high Curie temperatures 
in some magnetic semiconductors with direct band gap, such as Gai-a;Mni:As. 

PACS numbers; 75.50.Pp, 72.80.Ey, 75.30Hx, 75.50Dci 



I. INTRODUCTION 

The advantage of ferromagnetic semiconductors (FS) 
as a source of spin-polarized carriers is that they can be 
easily integrated into semiconductor devices When 
discovered, ferromagnetism at room temperatures with 
full polarization of itinerant carriers will be a major 
breakthrough in semiconductor electronics. Most the- 
oretical and experimental efforts have been concentrated 
on III-V, group IV, and Il-VI-based diluted magnetic 
semiconductors (DMS). These semiconductors are alloys 
in which some atoms are randomly replaced by magnetic 
atoms, such as Mn^+. 

Ferromagnetism in diluted magnetic semiconductors 
(DMS) is thought to be well understood in terms of the 
so-called p-d exchange model, which was first considered 
over 50 years ago 0,0,01 (See also Ref. 0,0,0 for a re- 
view). At concentrations of impurities above the Mott 
limit, i.e., as soon as carriers become delocalized, con- 
ventional model of FS is fairly simple. According to, for 
example, Ref. , the interaction between charge carriers 
in a semiconductor and spin-S impurities can be written 
as: 

U ^ - j s{r)Y,S,.P\r-Ri)d'r (1) 

i 

= -jpto^y^^ Sj, 

i 

where 5'^ and Ri are the spin and the position of an i- 
th atom of magnetic impurity, J^'^q = / J{r)(Pr, and 
we have used the fact that magnetic impurities are ran- 
domly distributed in the sample. A simple analysis 
then shows that, in the presence of one type of charge 
carriers, Curie temperature, Tc, is proportional to con- 
centration of magnetic impurities, Ni, and the square of 
the strength of the exchange interaction, a: 



n,5(5 + l)(jPlo)^Xo 



(2) 



Here fj,o is the magnetic moment of charge carriers, xo is 
the Pauli term in the spin susceptibility in the absence 
of impurities. At sufficiently large carrier densities, the 
spin polarization of the charge carriers at T = is given 
by: 



(3) 



In general, for any carrier density, the spin polarization is 
given by Zeeman-split Fermi surface, with the difference 
in chemical potentials for "up" and "down" spins given 
by: 



MT - Mi 



(4) 



The spin polarization at any filling is then very easily 
calculated from this equation. For example when the 
carrier density 



Tie < Tie 



{2JlUn,Sm*fl^ 
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(5) 



the carriers will be fully polarized at T = 0. When rie > 
ric, they are no longer fully polarized. The polarization 
of carriers is then determined by a parametric equation 
(with /i as a parameter): 



(net - n^i), 



where 



(2m> ± Jli^n,Sm*fl^ 



67r2 



(6) 



(7) 



and rie = riei + rie-f- 

Despite the simplicity of the basic concept of Zener fer- 
romagnetism, calculations of Tc for real materials become 
rather involved, and depend crucially on details of the 
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band structure, the p — d exchange matrix, or direct anti- 
ferromagnetic exchange between Mn^+ spins. Most the- 
oretical and experimental efforts on magnetic semicon- 
ductors have been concentrated on finding a new DMS- 
based material that would have a ferromagnetic transi- 
tion above room temperature, and would be possible to 
incorporate in thin film form with mainstream semicon- 
ductor device materials. There are theoretical predic- 
tions for Tc's above room temperatures in several classes 
of these materials |9llTol|. However, experiments indicate 
that the growth of Curie temperature with concentration 
of magnetic impurities saturates at 5-10% Mn doping in 
most ferromagnetic semiconductors; it may even start to 
decrease at higher concentrations of Mn. Since, accord- 
ing to the mean field p-d model, the Curie temperature 
Eq.@ grows linearly or faster (if the carrier concentra- 
tion changes) with the growth of Mn concentration, it is 
important to understand what limits the growth of Tc, 
and how this can be avoided. 

Experimental effort in FS concentrated mostly on 
Ini_^Mn:rAs'lll (Tc - 35K) and Gai_:rMn:rAs[ll [H 
0,^] {Tc ~ llOi^). While room temperature ferromag- 
netism in FS remains a theoretical possibility, the highest 
Curie temperature reported in the homogeneous sample 
with Mn concentration around 5% is 116K for Mn^-Gci^^ 
|l6j |. The saturation Tc growth with increased Mn con- 
centration is usually ascribed to increased disorder. Dis- 
order effectively introduces an exponential cutoff for the 
RKKY interaction and reduces its range 12). However, as 
we discuss in section VHI, weak disorder docs not change 
Curie temperature. Another possible reason is that the 
strength of direct antiferromagnetic exchange grows as 
the average distance between impurities becomes shorter, 
which, in turn, lowers the Curie temperature. Room tem- 
perature ferromagnetism has been observed in a number 
of compounds with large concentration of Mn impurities 
(see Ref. Q for a review), or, due to phase separation 
and formation of nanoclusters of these compounds ( such 
as MniiGeg in Ref.0|). However, large carrier concen- 
tration in these compounds limits the degree of spin po- 
larization that is necessary for device applications. 

Some theoretical studies ^3 claimed that the reason 
Curie temperature can't get higher above certain concen- 
tration of Mn is that the p—d model cannot be treated in 
the mean field approxination. However, since the RKKY 
interaction has a large radius, mean field treatment in 
absence of strong disorder should be justified 18]. 

In what follows, we adopt the idea that the reason 
for the discrepancy between mean field theory and ex- 
periment lies in the presence of additional interactions 
in the effective spin Hamiltonian of the problem, which 
are not present in the conventional formulation of the 
p—d model. To derive effective spin Hamiltonian, includ- 
ing additional interactions, we start from a more general 
Anderson model of deep magnetic impurities in FS. We 
show that, in case of a deep Anderson impurity in a semi- 
conductor, an additional long-range indirect exchange in- 
teraction appears in the effective Hamiltonian, which, if 



antiferromagnetic, severely limits Tc-s in these materials. 
On the other hand, for direct gap semiconductors, this 
interaction could change sign and become ferromagnetic, 
if the magnetic impurity lies deep enough. We derive 
conditions under which this becomes possible. As a re- 
sult, ferromagnetic correlations would become enhanced, 
not reduced, as the concentration of magnetic impuri- 
ties grows, and high-temperature ferromagnetism could 
be possible even without carriers. The indirect exchange 
between two deep impurities, whether ferro- or antifer- 
romagnetic, is stronger than the RKKY interaction, and 
thus could produce large Curie temperatures. For ex- 
ample, it could provide a possible explanation of high- 
temperature weak ferromagnetismfioj in LaajCai^ajBg or 
recently discovered CaB2C2'20] (CaB2C2 has Tc = 770K 
and M = lO^'^fis)- Experiment 20.. .2Lj indicates that 
these materails are direct band semiconductors with a rel- 
atively small band gap, and that impurities play a major 
role in establishing the new high-temperature ferromag- 
netic state. However, other reportsj23| claim that high- 
temperature ferromagnetism in these materials is not a 
bulk effect. Rather, it is related to clustering or new 
boron phases with Fe or Ni magnetic impurities. Under- 
standing the mechanism of weak ferromagnetism in these 
materials could lead to a discovery of more members of 
this class. In either case, as we report in this paper, there 
is a theoretical possibility of a new type of FS with direct 
band gap, where ferromagnetism is not carrier-driven. 

This paper is organized as follows. In Section II we 
discuss the general procedure of the derivation of the ef- 
fective low-energy Hamiltonian for Anderson impurities. 
In section HI we discover that, due to strong hybridiza- 
tion, an Anderson impurity is no longer a local center; 
its localized wave function acquires a finite range, which 
is directly related to the corresponding "Bohr" radius of 
a charged impurity. Sections IV and V are devoted to 
the derivation of the effective exchange Hamiltonian for 
magnetic impurities. In Section VI we explore the conse- 
quences of the large range of interaction for magnetism, 
such as a high Curie temperature in case of a direct band 
gap if magnetic impurities are dense enough. In sec- 
tion VII we consider a minor modification of the effective 
Hamiltonian in case of higher spin (such as 5' = 5/2 for 
Mn), and the application of these ideas to Gai^j^MuajAs. 
In section VIII the influence of weak disorder and interac- 
tions on Curie temperature is briefly discussed. Section 
XI provides a summary and conclusions. 



II. THE EFFECTIVE HAMILTONIAN 

We start by considering a simple model of FS. In III- 
V systems, such as Gai_a;Mna;As, it is well established 
that the Mn ions substitute for Ga, and contribute itiner- 
ant holes to the GaAs valence band. Experimentally, the 
hole density is typically a small fraction (15% or so) of 
the Mn concentration, perhaps due to strong localization 
of carriers on Mn and other defects, so Gai_xMn2;As can 
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be considered partially compensated. The Mn ion has 
half-filled d-shell, which acts like a spin-5/2 local mo- 
ment. The Anderson model, which is more general than 
the p — d model, should completely account for all the 
physics of FS. It is well known that, when spins are well 
localized, the single-impurity Anderson Hamiltonian is 
reduced to the p — d Hamiltonian by the Schrieffer- Wolff 
transformation ^23j . For many impurities, this may no 
longer be the case. Let us start by considering a single- 
orbital Anderson Hamiltonian: 



where 
Ho 



'tpcr 



(8) 
(9) 



E 



^dl^^dna + Ud^^dl I djii dii'\ 



Here the first sum (p) is taken over the reciprocal 
space, the second sum (n) over real space impurity sites. 
U is the on-site Coulomb repulsion term. Typically, U 
is very large (~ 5eV), and can be taken to be infinite. 
Here (p) are the energy band spectra of conduction and 
valence bands {i = 1, 2). 

The hybridization term in the model Hamiltonian, Hy, 
accounts for the p — d hybridization between impurity 
sites and conduction and valence bands: 



7V1/2 



ip-Rn 



h.c 



} (10) 



pnai 



This model is a reasonable generalization of the p — d 
exchange Hamiltonian, usually considered in the litera- 
ture. Because of large on-site Coulomb repulsion, the 
d-levels are half-filled. While we consider the case of a 
single d-orbital, a generalization to S = 5/2 Mn ion is 
straightforward (see section VII below). 

The Anderson Hamiltonian Eq.© describes a very 
complicated problem. However, under the U — oo con- 
straint, it can be reduced to the problem of Heisenberg 
spins (in case of a single d-orbital, spin-1/2). The low- 
energy effective Hamiltonian is equivalent to Eq.(|Sll in 



the limit ksT <C A^, where is the energy difference 
between the impurity d-level and the top of the valence 
band, or the energy difference between the impurity d- 
level and the bottom of conduction band. Since typical 
gap values in semiconductors arc of the order of leV, this 
is usually a valid assum ptio n. 

Following Refs. 0123, the effective spin Hamilto- 
nian can be derived by expanding the S-matrix (or, at fi- 
nite temperatures, the partition function) in V-s, and re- 
expressing various time-ordered processes in terms of spin 
operators. Then, these processes are collected back un- 
der exponent, to obtain the effective Hamiltonian. This 




FIG. 1; The second-order contribution to the effective Hamil- 
tonian, shown as a time-ordered process. 



method allows one to obtain consistently the interaction 
between spins and carriers (electrons or holes) in conduc- 
tion and valence bands, and carrier-carrier interaction. 
In a way, the concept is similar to perturbative renor- 
malization group, since we arrive at a low energy effec- 
tive Hamiltonian by integrating out higher-energy states. 
Treating hybridization term in Eq. (|SJ) as a perturbation, 
we can rewrite the partition function as 



Z = Tr[exp{-(3Ho)S{P)], 



(11) 



where 



J 



Si(3)=Texp l-j^^Hv{T)dT \ = Texp(- ^ drF, {vI/t^,(T)d„,(T) + dL(T)vE'„,,(T)} 



(12) 



The problem of finding the effective Hamiltonian is then to reduce these expressions to the form: 



Z = rr(exp[-/3He//]), 



(13) 
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using the fc^T <C condition, i.e., the fact that the 
local levels are almost always occupied, and that transi- 
tions to conduction and valence bands are absent at low 
temperatures. Various terms in the effective Hamiltonian 



The order of operators dl^^{T2)dna' {ti) (t2 > ti) is fixed 
by the assumption of strong Coulomb repulsion on the 
n-th center. Note that, in the second order, because all 
centers are filled, only one band contributes to the effec- 
tive Hamiltonian. The "filled" band does not contribute, 
because of the Pauli principle (this is not the case at a 



The first term in Ea. H15(l is a c- number. In the limit 
Ai ^ ksT it is possible to put ti ~ T2 in the second 



where Sn is the localized spin of the n-th impurity. The 
first term corresponds to the renormalization of the en- 
ergy of the localized level, while the second term involv- 
ing spins of localized impurity and carriers is nothing but 
the ordinary Hpd, the p-d model discussed in section I. 

The next order terms in the effective Hamiltonian are 
fourth order in V-s. There are sums over two local cen- 
ters, m and n in S'4(/3). Also, there are contributions 
to S4{(3), which we will denote S'^{/3), that are already 
accounted for in the effective Hamiltonian Ea. l|16|l : 

S'M = \ j j dridT2T{Hpd{Ti)Hpd{r2)}. (17) 

These contributions need to be subtracted from 6*4 (/3), to 
get the fourth-order (in Vi) contributions in the effective 
Hamiltonian. The most important contribution is the 



can then be associated with certain time-ordered virtual 
processes. For example, the first non-zero contribution 
to Heff is from the second order term in the expansion 
of Eq.linil in V,: 



(14) 



finite U). The perturbative time-ordered process corre- 
sponding to the term in Ea. H14|) is shown in Fig^ 

We may rewrite VPriCTi (t'2)^1o-'i ('''1) interaction 
representation as: 



(15) 



term in Eg. 1)15(1 . Thus, for the second-order contribution 
in the effective Hamiltonian we may write: 



(16) 

I 

effective exchange interaction between localized spins: 



H^'^'^ — — ^ J{Rn — Rm)S„Sm (18) 

n^rn 

This exchange interaction is the result of 2 time-ordered 
processes shown in Figs. |2 121 One is the superexchange 
(Fig. 121) , which is a result of the localized spins exchanged 
through the empty conduction band. The other pro- 
cess is the Bloembergen- Rowland termpGf (Fig. EJ, an 
exchange process through both conduction and valence 
bands. The form of these contributions will be discussed 
in greater detail in sections IV and V below. 

In addition, other interesting contributions arise as a 
result of fourth-order processes, such as p-d scattering by 
spins on two centers. 



Jo Jt, 



I 

*«.i(t2)4'L'i(^i) = {*nal(^2)*L'l(^l)}+ - *L'l(n)*"-l(^2). 

I 



I 

p ^ naa' ^ 



J 



^SL = X2 E *(^" - Rrn)'^la,,i^nal[SnS^5,,. + Z ( [5„5„]<3-,.,)] , (19) 

^1 -I. I 
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nontrivial local contribution, 



(4) 
local 



A3 



E 



* nal ^ np'l 



■ ^p' p + {o'p' pSn 



(20) 




FIG. 2: The fourth-order antiferromagnetic superexchange 
contribution to the effective exchange interaction between two 
localized impurities. 



and corrections to the energy of the local level and 
ground state energy, which are dropped. These nontrivial 
terms, however, are higher order in carrier density, which 
is small in magnetic semiconductors. Thus, they don't 
play any significant role in magnetism, and, hence, will 
not be discussed in any detail. We refer the reader to 
references 24^ 25j 2^ : where these terms were discussed 
in detail in connection with the 3-band model of cuprate 
superconductors. 



III. RENORMALIZATION OF LOCAL LEVEL 
DUE TO HYBRIDIZATION. 

To clarify the new physics of an Anderson impurity, 
in comparison with p — d magnetic impurity, let us first 
consider a single-level problem (i?„ = 0), interacting only 
with the conduction band (i.e., we take V2 = for sim- 
plicity). The Hamiltonian is defined on the manifold of 
wavefunctions {$0, ^ipu ^lp2^ which correspond to 
the exact solution for a center in crystalline lattice in 
the absence of hybridization. These wavefunctions form 
a complete orthogonal basis for the single-particle prob- 
lem. We may now introduce the hybridization as a per- 
turbation into the Schrodinger equation. In the matrix 
representation, the single-level Hamiltonian takes a very 




FIG. 3: The fourth-order Bloembergen-Rowland term in the 
effective exchange interaction. 



simple form: 



eo vi vi 
vi eipj 
vi eip2 



(21) 



Here vi = Vi/N^^"^. This Hamiltonian can be easily 
diagonalized, by looking for a solution for the ground 
state wave function as a linear combination of all single- 
particle states: 



(22) 



where the coefficients obey the following set of equations: 

{^o-e)xo + viJ2xp = (23) 
p 

ViXo + (eip - e)Xp = 

Solving Eq. (|23|l gives the new position of the localized 
level: 



e = eo 



(24) 




FIG. 4: Localized centers acquire a finite radius in the Ander- 
son model. This radius can be as large as 2-3 lattice spaces in 
magnetic semiconductors. An overlap of wavefunctions from 
two localized level leads to the new physics, which is not cap- 
tured by the traditional p — d model. 

Thus, in presence of hybridization the position of the im- 
purity level is shifted. If the depth of the level Ai Di, 
where Di is the bandwidth of the conduction band 
(which is usually a valid assumption for semiconductors, 
where the band gap is much smaller than the correspond- 
ing bandwidths of conduction and valence bands), the 
sum in Ea. (|24|l can be performed numerically. The cor- 
rection to the energy level, 

V? 

e-eo--^, (25) 

turns out to be small for small enough hybridization pa- 
rameter Vi ^ \/DiAi. It is responsible for a slight down- 
ward shift of the energy level. Such corrections will thus 
be dropped in our further discussion. 

The wavefunction of the localized level, according to 
Ea. (|23|l . acquires an admixture to its decay in the form: 

r 

Here a is the lattice constant (we assume a simple cubic 
lattice), TOi is the effective mass of carriers in the empty 
band (for hole-doped materials, the valence band is the 
"empty" band, since it is empty of holes), Ai is the en- 
ergy difference between impurity levels and the bottom 
of the empty band. Thus, a new length scale enters the 
problem, 

Ro = (2miAi)-i ~ a^Di/Ai, (27) 

which, while small compared to the average distance 
between doped carriers, may considerably exceed inter- 
atomic distances (the large parameter is the square root 
of the ratio of the bandwidth to the energy gap in a 
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FIG. 5: Positions of conduction and valence bands. 

semiconductor), and become comparable to the average 
distance between magnetic impurities. The overlap of the 
localized wave functions gives rise to a new contribution 
to the exchange integral, which will be derived below. 
The picture of finite-radius localized magnetic impurities 
with overlapping wavefunctions is shown in Fig0| 

Existence of a large (compared to lattice spacing) ra- 
dius for deep impurities was first noticed by Keldvshp^. 
who analyzed deep charged impurities in ordinary semi- 
conductors. In his case, however, this scale corresponded 
to ordinary Bohr radius for a deep impurity ( modified 
by the fact that in that case one had to consider Dirac 
Hamiltonian for the k — p model of semiconductors) . In 
our case, the large length scale comes from the hybridiza- 
tion of the localized impurity level with conduction and 
valence bands. We will also see below in Section VI that 
effective length scale, at which direct interaction between 
impurities starts to matter, grows logarithmically with 
decreasing temperature in the framework of percolation 
theory or virial expansion. Thus, at low temperatures, 
such as Curie temperature, the effective range of interac- 
tion between two impurities becomes even larger. 

IV. THE EXCHANGE HAMILTONIAN. 

The problem of magnetism in the Anderson model can 
be reduced to the spin Hamiltonian Eq. (|18|l , described by 
the processes shown in FigsHJO In case when the empty 
band (the "empty" band can also be the valence band, 
in case of hole doping) gets partially filled by carriers 
(electrons or holes) , the time-ordered processes shown in 
Figs m include all possible contributions - the familiar 
RKKY exchange interaction, the superexchange, and the 
Bloembergen-Rowland interaction. In presence of carri- 
ers, the superexchange gets modified by the RKKY ex- 
change term. We shall see below, that the Bloembergen- 
Rowland interaction also gets modified in presence of 
dopants. This modification is not accounted for by the 
RKKY interaction. As we have discussed in the previous 
section, localized impurities in the Anderson model ac- 
quire a large radius (compared with the lattice constant, 
a). Thus, exchange interaction between localized spins 
is, in general, long range. This would allow one to use 
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mean field to describe the ferromagnetic ordering, if the 
concentration of localized spins is high enough. On the 
other hand, if the concentration of magnetic impurities is 
low, one can use percolation theory or virial expansion. 

Let us first consider exchange integrals given by the 
processes shown in Figs El El for a semiconductor with 
conduction and valence bands separated by some general 
reciprocal space wavevector Q, as shown in FigEl When 

I 



J(q) 




FIG. 6: Superexchange contribution to exchange integral in q- 
space; modification by RKKY in presence of carriers is shown 
by a dotted line. 



Here i, j are the band indices for conduction and va- 
lence bands, ji is the chemical potential, n(e — /x) is the 
Fermi-Dirac distribution function. The empty band cor- 
responds to i — 1, while the filled band has the band 
index i — 2. We can see now, that for T <^ Aq, where 
Aq is the band gap, the contributions J21 and J22 are 
absent. Jn is the superexchange, which corresponds to 
the process shown in Fig. El In presence of carriers, this 
process, and the corresponding expression, also includes 
the RKKY contribution. When carriers are present, the 
effective Hamiltonian can be written in terms of impu- 
rity spins only, with carriers " integrated out" . In this 
case, impurity spins interact with the exchange Hamil- 
tonian Eas. (P5|) . lf^ . which includes RKKY. If, however, 
we wish to retain carriers, and the carrier-impurity p-d 
exchange interaction, as we described in section II, then 



Q = 0, it is a direct band gap semiconductor. Surpris- 
ingly, all exchange processes can be written in a relatively 
compact way: 

J{R)^J2J^^(R)' (28) 

where 



(29) 

I 

J(q)' 




FIG. 7: Bloembergen-Rowland contribution to exchange in- 
tegral in q-space; modification in presence of carriers is shown 
by a dotted line. 



the RKKY part wiU be the part S'{P) in the 4-th order 
(see Ea. ((T7|l '). which has to be subtracted, and the indi- 
rect exchange interaction will be described by Eq. (|29|l at 
zero doping, i.e., with the factor [1 — n(epi — /i)] in Ea. l|29|) 
replaced by 1 . These two descriptions are only equivalent 
when the chemical potential for doped carriers, /x <C Ai. 
Basically, when the doped carriers are retained, we "in- 
tegrate out" electrons and holes in the original Anderson 
Hamiltonian up to the scale /x. When we consider the 
Hamiltonian for localized spins only, we integrate out all 
carriers. It makes sense to define the RKKY interac- 
tion for the Anderson model as the difference between 
the doped and undoped cases in Eas. lpSll . lt^ . When 
fi ~ Ai, the Anderson model cannot be reduced to the 
p — d model, except near the Fermi surface. The Ander- 
son model expression for the RKKY interaction in real 



' (27r)6 



(fpcfq 



[1 - njep, - /i)]e'q« 
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space will then differ from the corresponding p — d model 
expression at short distances, which are important for 
ferromagnetism. We won't consider this subtle point any 
further, since we will always assume the carrier density 
in FS to be small, i.e., /i ^ Ai. In this section and the 
next section we will deal with an effective Hamiltonian for 
spins only, i.e., when all carriers are integrated out. This 
procedure for the Anderson model is always justified (al- 
though, as we commented above, a reduction to the p—d 
model is not). J12 is the Bloembergen- Rowland interac- 
tion, which also gets somewhat modified by doped car- 
riers. Note that the expression for the exchange integral 
Ea. (|29|l . obtained from the Anderson model, differs from 
the original expression of Bloembergen and Rowland ; 
it includes the term (epi — eo)^ in the denominator, which 
makes the exchange integral in the Anderson model long 
range. Both contributions to exchange integral can be 



J 



ViO^mipF 



evaluated numerically and analytically; they are shown 
schematically in Figs and [7| 

We see that the superexchange is, in general, antiferro- 
magnetic, while the Bloembergen-Rowland contribution 
favors a spin density wave with a wavevector Q, sep- 
arating conduction and valence bands. Both exchange 
processes are, in general, of the same order, and their 
relative strength is determined by the corresponding hy- 
bridization parameters Vi and V2. The superexchange 
integral in the q - space, in the absence of carriers, can 
be easily found analytically: 



.hi{q) = - 



1 



7rV2miAi + 8miAi 



(30) 



In presence of carriers, Jn {q) gets modifed by the RKKY 
interaction: 




(2^)2A? 



(31) 



where Pf = ^/2miJI is the Fermi momentum for doped 
carriers. The superexchange and RKKY contributions 
can be easily rewritten in real space: 



and 



with 

Jse{R) 



Jii{R) — Jse{R) + Jrkky{R), 



(32) 



exp{-V8TOiAii?}, (33) 



and a standard expression for the RKKY interaction, tak- 

Jpd 



ing into account that Jpd = 2F]^/Ai (note that in Ea. l|18() 



we sum over each pair of impurities twice) : 

Jrkky{R) = -JrkkyF{2pfR), (34) 

where 

, cos(a;) sin a; 
Fix) = - — ^, (35) 



J, 



RKKY — 



Via°mipp 



(36) 



We define the RKKY interaction in the Anderson model 
formulation as the difference between Jii(i?) in doped 
and undoped cases, given by Ea. H29(l . As we empha- 
sized above, since it depends explicitly on impurity en- 
ergy level, when Ai ~ fi, this expression differs from 
the usual RKKY form. Its asymptotic behavior at large 
distances, certainly, docs not change. In what follows 
we assume that the carrier concentration is always low, 
H <^ Ai, so that Jrkky{R) is given by the standard 
expression. Note that if Jpd and the level position are 
known, the exchange integral Ea. (|32|l contains no free 
parameters: 



Jn{R) = -Jrkky F{2pfR) 



3/2 



7r(TOi Ai) 
2^p'fR 



,-2\/2miAi_R 



(37) 



Similar to superexchange, the Bloembergen-Rowland in- 
teraction can also be written in two parts. The first part 
is from the empty and filled bands, while the second part 



takes doped carriers into account: 

•/12(g) = JsRiq) + ^(2(9)- 



(38) 
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The first contribution in Ea. (|38|) favors a spin density 
wave witli a wavevector Q, separating conduction and 



valence bands, and has a rather cumbersome form: 



JBR{q) 



a?ViV2m\'m2 
7r\/2miAi 



miAi + maAi + |- + maAo - ^2(mi + m2)Ai J2m2Ao + g^-i 



l+m2 



Af (mi + 7712)^ + Ai(q2(rni - 7712) - 2to2(™i + m2)Ao) + + TO2Ao^ 



r 



(39) 



Here and below in this section the notation will be 
slightly different; q is now the difference from Q (the 
wavevector separating the bottoms of the two bands), 
i.e., we change our notation in the following way: q — 
Q — > q. In general, the Bloembergen-Rowland contri- 
bution Ea. H39|l is quite cumbersome, which makes it im- 



possible to find an explicit analytical expression in real 
space, except in some limiting cases. However, the lead- 
ing asymptotic behavior of Jbr{R) at large distances can 
be derived. If the deep impurity level lies inside the gap, 
Ai < Ag, we find an exponential decay of exchange cor- 
relations: 



J 



Jbr{R) 



a^ViV2mim2 




cos {Q ■ R)e-^/^' 



(40) 



where the range of the integral is: 

1 

Ri = 



V2miAi + V2to2(Ao - Ai) ' 



(41) 



On the other hand, when the deep impurity level enters 
the filled band, Ai > Aq, the Bloembergen-Rowland con- 
tribution has an oscillating decaying asymptotic: 



Jbr{R) — 



a^ViV2mim2 



mi R 
cos ■ 



2Ai 



m2 



R 



2(Ai 



sin— ) coslQ ■ R)e-^/^-^° 
Ao Rci ' 



(42) 



r 



where 

Rci = 



^2mi(Ai- Ao) 



, R, 



cO 



1 



\/2toiAi 



(43) 



The contribution from doped carriers to J12 can be 
derived as well: 



^{2(9) 



In real space, this contribution becomes: 



(44) 



We see that if Q = 0, i.e., when our FS has a direct 
band gap, the Bloembergen-Rowland mechanism gives a 
large ferromagnetic short-range contribution to the ex- 
change integral. At high enough concentration of impu- 
rities, ferromagnetic properties are then detcrmimed by 
the value of exchange integral at q = 0. This value can 
be easily written down (it is only valid for a direct band 
gap FS): 



JM - COS (Q . i^)e-^»^« (45) 



J(q = 0) = Jn(q = 0) + Ji2(g = 0) 



V^m2 



7rV2miAi \ 8A1 (^ mi -|- m2)Ai 4- \/2m2Ao)^ 



(46) 
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In this section we have derived exphcitly the interac- 
tion between two magnetic impurities in a semiconductor 
with one conduction and one valence band. Typically, the 
band structure of semiconductors is more complex than 
that (for example, GaAs has light and heavy hole bands). 
This situation is considered in the next section. 



V. WHAT IF A SEMICONDUCTOR HAS MORE 
THAN TWO BANDS? 

When a semiconductor has many bands, the analysis is 
just as straightforward as in case of two bands considered 
in the previous section. The exchange integral between 
two impurities is still given 

by Eqs||2Hll,(I2ni), where the 
sum now goes over all pairs of band indices. The case 
when i = j is an "empty" band index corresponds to 
a superexchange contribution (Fig|2J), considered in de- 
tail in the previous section; — when z is a "filled" 
band index, because of the Pauli principle; Jij is the 
Bloembergen- Rowland contribution (FiglJJ, when i is an 
"empty" band index, and j is a "filled" band index. The 
contribution of a new type appears when one has two 
or more "empty" bands in a semiconductor. Then there 
will be a superexchange contribution of the type shown in 
FigI21 where now the carriers are exchanged through two 
different "empty" bands - the process (1) in Fig|5|puts a 
carrier from one magnetic impurity into the first "empty" 
band, while the process (2) puts a carrier from the other 
magnetic impurity into the second "empty" band. Natu- 
rally, there are two such contributions, Jij and Jji, where 
i j indices correspond to two different "empty" bands. 
The exchange Hamiltonian in any FS is a sum of all pair- 
wise contributions listed above. 

Let us now consider the new two-band superexchange 
contribution in detail. In the most general case, the bot- 
tom of the second "empty" band is shifted from the bot- 
tom of the first "empty" band by a wavevector Q. The 
bottoms of the two bands also lie at two different posi- 
tive energies Ai and A2 relative to impurity level, unless 
there is a symmetry-related degeneracy, and the carri- 
ers in the first and second bands have different effective 
masses mi and m2. 

Let us first consider the simplified case when the bands 
1 and 2 are identical (i.e., Ai = A2 and mi = ^,2), but 
separated by a wavevector Q. Then we don't have to do 
the calculation, since ep2 ~ ep_Qi, and it can be easily 
seen from Ea. lf^ . that 

Jl2(q) = Jll(q-Q), (47) 



and 

Ji2(i?) = Jii(i?)e'«^. (48) 

Summing all contributions from 2 identical empty bands, 
since Jii(i?) — J22{R), we get: 

J2bse = 2 Jii (i?) [1 + cos (g -R)]. (49) 
This also includes RKKY-type contribution for two iden- 
tical bands at finite doping, described by: 

J12RKKY + J21RKKY — "^JRKKviR) cos {Q ■ R), (50) 

and the ordinary RKKY interaction from both bands. It 
is rather obvious from Ea. H50|) . that at a finite density 
of carriers there are ordinary one-band RKKY contribu- 
tions from both bands. The two-band RKKY contribu- 
tion oscilates much more rapidly in space, as cos(<5i?)), 
unless Q — 0. In the latter case, the two-band RKKY 
contribution just enhances the contributions from the two 
separate bands. In the case when Q = 0, the 12 and 21 
exchange integrals take the same form as shown in Fig|Sl 
For Q 0, these contributions have the same form in 
q-space as for Q = 0, centered at the wavevector q ^ Q. 

Now that the form of new two-band contributions to 
the exchange integral has become clear from the sim- 
plified case, let us consider a more general case of two 
completely different "empty" bands. The bottom of the 
second band is shifted from the bottom of the first band 
by the wavevector Q. Of course, 12 and 21 contributions 
to J{q) will still be of the form shown in FigEI with the 
bottom dX q — Q, but the corresponding expressions be- 
come more combersome. Let us assume that Ai < A2, 
i.e. the bottom of the first band lies below the bottom of 
the second band, and the first band gets filled by carriers 
first. We consider three separate cases, (1) both bands 
are empty; (2) the first band gets filled by carriers, but 
the second band is empty, and (3) both bands get par- 
tially filled by carriers. When Ai = A2, i.e., the two 
bands are symmetry-related, we only have cases (1) and 
(3). As before, we assume that the carrier concentration 
is very small, i.e., the Fermi energy for doped carriers is 
much smaller than any other energy scale in the problem, 
except temperature. 

(1) Both bands are empty. 

Let us now consider the exchange interaction arising 
from two empty bands. The integral in Ea. (|29|) can be 
easily calculated: 



TO2 



7i2+2i(9) - - _ Q)2 + (^2^1 Ai + V2m2A2)^ Wa7+VA^ ^ ^ 



I 

In the coordinate space, it takes the following form: where the range of this interaction is given by: 

Ji2+2i(H) = -V^i2Cos(Q.H)iexp(--^), (52) ^12= /t^^I /^^^ ' (^3) 

R R12 V2miAi + v2m2A2 
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and 



Vl2 = 



^1 ^2 " '"'l"^2 



(54) 



(2) One hand gets filled. 



Let us assume in this section that band one gets filled 
first, i.e., Ai < A2. We also assume that the chem- 
ical potential (counted from the bottom of band one), 
/ii <C A2 — Ai, i.e., the number of carriers is low. Of 
course, the main part of the exchange inegral is still 
given by Ea. (|52|l . However, at finite doping there are 
corrections. Obviously, since only band one gets carriers, 
(5J21 = 0. At small filling, RKKY-like correction can be 
easily calculated: 



5Ji2{q) 



3^2A2(2^2[A2- Ai]+g2 



(55) 



Here q denotes the difference from Q, the wavevector 
separating the bottoms of the two bands. In coordinate 
space it can be written as: 

6J,2{R) - ^'^ffip -sQ . H, (56) 



where 



Ra 



V2m2(A2- Ai) 



(57) 



(3) Both hands get filled. 

When Ai = A2 = A by symmetry, both bands get 
filled by carriers simultaneously. The result of this is 
an RKKY-like correction to exchange integral Ea. H52|l . 
which is given by Eq. (|50|) for two identical bands. If the 
bands have different masses, the interaction is RKKY- 
like long range, but the general expression is quite messy. 
Here we give the its value at g = Q, which is important 
for ferromagnetism, when Q = 0: 



SJi2+2iiq = Q) 



2V{V{a-^^/2jlmim2 
7r2A2(^ 



(58) 



To summarize this section, the exchange interaction 
in any magnetic semiconductor is composed of pairwise 
contributions: the superexchange contribution through 
one "empty" band, the Bloembergen-Rowland contribu- 
tion from one "empty" and one filled band, and a pair- 
wise contribution from two different "empty" bands. The 
last contribution was considered in this section. It in- 
cludes two-band long-range RKKY-like interaction at fi- 
nite doping, which is usually not taken into account in 
the p — d Hamiltonian. Note that when Q ^ 0, the two- 
band RKKY interaction will favor spin glass order, not 
ferromagnetism, since it oscillates very rapidly in real 
space. 



VI. THE MEAN FIELD APPROXIMATION, 
VIRIAL EXPANSION, AND PERCOLATION. 



Now that we have obtained all terms in the exchange 
interaction Ea. (|29|l . we can proceed to calculate mag- 
netic properties of ferromagnetic semiconductors. In 
general, the RKKY exchange interaction is always long- 
range. The shorter-range contributions may or may not 
be treated in mean field, depending on the ratio of their 
range to the average inter-impurity distance. We as- 
sume here, for simplicity, that one of indirect exchange 
processes considered in the previous two sections domi- 
nates the short-range physics. In case of FS with indirect 
band gap, the Bloembergen-Rowland exchange contribu- 
tion favors antiferromagnetism with a lattice wavevector 
Q. It oscillates very rapidly in real space. Since im- 
purities are distributed randomly, it would favor a spin 
glass ordering. The contributions at Q = 0, which influ- 
ence ferromagnetism, are the RKKY and the superex- 
change. From Ea. H37() one can easily see that, when 
impurity concentration rii ^ (rriiAi)^^^ , the superex- 
change contribution totally suppresses ferromagnetism 
brought about by the RKKY exchange interaction. On 
the other hand, when Ui <^ (miAi)'^/^^ the superex- 
change at an average distance between impurities is sup- 
pressed, and the RKKY interaction gives rise to ferro- 
magnetism, as in the ordinary pd-model. Corrections 
due to superexchange can then be calculated using the 
virial expansion approach j29|. We can rewrite the domi- 
nant exchange contribution from Eq. (|37|) in the following 
form: 



J{R) = -JrkkyF{2pfR) 



R/Ro 



(59) 



where Vq = t^Jrkky l{&'^P%Rt)^ Rq = l/V8miAi, in 
case when it is given by one empty band only. If contribu- 
tions from more than one "empty" band are important, 
Vq and Jrkky are, in general, unrelated to each other. 
The reason for this is that the largest short range con- 
tribution typically comes from the lightest bands, which 
would produce the longest range indirect exchange (if all 
hybridization parameters are of the same order) . On the 
other hand, the dominant RKKY contribution is the one 
from the heaviest band, since the carriers in that band 
have the largest Fermi wavevector; there will also be ad- 
ditional RKKY contributions from many bands, which 
were analyzed in the previous section. In what follows 
we consider the simplest case, when only one heavy band 
is relevant for RKKY. 

In case of a FS with a direct band gap, the 
Bloembergen-Rowland mechanism gives rise to ferromag- 
netism, and we can either have a ferromagnetic or an 
antiferromagnetic short-range exchange contribution, de- 
pending on the relative strength of corresponding ex- 
change processes. If the Boembergen-Rowland term 
dominates the physics at short distances, the impurity 
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spin Hamiltonian becomes: 



the number of carriers: 



J(i?) = -JrkkyF{2pfR) + Vbr^ e-^/«^ , (60) 

In general, the Bloembergen- Rowland term Eq.(jSHIl also 
has an antiferromagnetic contribution from doped carri- 
ers. It then can be rewritten in the following form: 



R/Ri 



Vdo 



Rg 

Rn 



-R/Rg 



(61) 



where Rg — l/-\/8m2Ao, Vbr. and Vdop are given by the 
coefficients in front of exponents in Ea. (|40() and Ea. (|45|l . 
However, since 



Vdop /Vbr 



(miAi)3/2 



« 1, 



(62) 



while Rg is still quite small, the contribution of doped 
carriers to Tc through the Bloembergen-Rowland mecha- 
nism can be neglected. Thus, if one exchange process 
dominates the physics are short distances, the short- 
range exchange contribution has the same form, and can 
only differ in sign. In what follows, we consider the gen- 
eral form of exchange integral Eq. 1)59(1 , assuming that the 
exchange constant Vq , its sign, and its range Rq are those 
of the dominant indirect exchange contribution. Note 
again that in our definition of exchange integral we sum 
over impurities twice. We will also consider the case when 
the short-range part in Ea. l59|) is rapidly oscillating, 

JsR ^-Vo^e-"^"" cos (Q-R). (63) 

As we have seen above, this happens when two bands 
separated by wavevector Q contribute the most to the 
short-range exchange interaction. Vq has a negative sign 
for the Bloembergen-Rowland contribution, and a pos- 
itive sign for the superexchange. When more than one 
process is important for short-range physics, the problem 
can always be solved numerically for a given set of pa- 
rameters. Here we obtain an analytical solution in several 
limiting cases. 

(1) Dilute system with almost no carriers. 

In the absence of carriers, the type of order and Tc is 
determined by the short-range part of interaction. If fer- 
romagnetic Bloembergen-Rowland interaction dominates 
at short distances. Curie temperature is approximately 
given by the ferromagnetic interaction taken at the aver- 
age distance between impurities @: 



Tc ~ 2.3VoS^Ron]^^ e 



Vo>0, 



(64) 



which is valid when rii <C V-Rq- A comparison with 
Eq.(|2Jl shows that this gives the following condition on 



PF = (STT^ne)^/^ < Pqf 
27.6Von^SRo 



(65) 



POF 



-0.87/(i?,oTii 



When the short-range interaction is antiferromagnetic, 
it favors spin glass order, and ferromagnetism is absent 
when there are no or almost no carriers. The condition 
on the number of carriers for ferromagnetism to be absent 
is then given by: 



PF < POSG 



27.6\Va\n^Ro ^_o.87/(i?on,'/') 



(66) 



Finally, for a rapidly oscillating short-range part Eq. I|63|) 
this condition can be rewritten as: 



PF < POSG 



27.6\Vo\n^R< 



to 



,-0.87/(flon-^^) 



For the Anderson Hamiltonian, J^^g ~ 2Via^/Ai. 
Note that if the ratio rie/ni is fixed, pp grows much slower 
with Ui than po- Thus, conditions in Eas. (|()^ . (|()()|l . (|(r7)l 
are likely to be satisfied at some finite concentration 
of magnetic impurities, ric <C <C ^/Rq- For exam- 
ple, for antiferromagnetic or oscillating interactions, Eqs 
(|65|l . (|67|) will define the concentation of impurities above 
which carrier-driven ferromagnetism disappears. 

(2) Dilute system with carriers. 

When the carrier concentration is large, i.e., pp 3> po 
in Eas. ((65)1 . H66 |l . ((67(1 . but the concentration of magnetic 
impurities is still small, rii <JC ^/Rq, short-range interac- 
tions between magnetic impurities will result in a correc- 
tion to Tc, which can be calculated by virial expansion 
(see, for example, Ref.p^l. Following Ref.(5ii], since the 
range of RKKY interaction is large, we can represent p—d 
interaction between carriers and magnetic impurities by 
a mean field Zener Hamiltonian: 



(67) 



H 



MF 



jpd 



(68) 



where s is the density of ordered spin of the carriers, 
which is assumed to be constant in space. In addition, 
there is a relatively short-range exchange interaction be- 
tween impurity spins, given by the processes described in 
Sections IV and V: 



He 



7 , J{Ri ~ Rj)SiSj 



(69) 



The short-range exchange integral does not include the 
carrier contribution, since it is already accounted for in 
the mean field Hamiltonian. In the Zener model Eq. 1(68(1 , 
Curie temperature Eq. ^ can be found by minimizing the 
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free energy density of the system of carriers and spinsj^, 



(70) 



with respect to s, and finding when the solution at small 
s first appears. Here 



2X0 ' 

and Fi is given by the usual Zecman term, 



(71) 



sinh [jP%s{S+ 1/2) /T] 

F, = -n,T In —-— (72) 

sinh[JlUs/{2T)] 

To find the virial correction to Eq.© from Ea. (|69|l . we 
need to include into F the contribution from two mag- 
netic impurities, when they are close enough: 

F2^ = \nf j d'R[F2^{R) - 2F,{R)], (73) 

and calculate it from the Hamiltonian 

Fl = HmF + Hexch- (74) 

This can be done, since the Zeeman Hamiltonian for two 
spins, interacting via direct exchange interaction, can 
easily be solved. The integral over this solution, how- 
ever, can only be taken with logarithmic accuracy. The 
contribution from two impurities is important when the 
distance between them is 



R< Rnln 



\Vo\/T 
ln(|K,|/T)- 



(75) 



Finding i<2i, and repeating the minimization over s, we 
obtain, as expected, that a for ferromagnetic exchange 
interaction (Vq > 0) Curie temperature is enhanced: 



ST^[F] AnS 

- 3(5+1) ln{\Vo\/T, 



(76) 



For an antiferromagnetic exchange interaction (Vb < 0), 
the Curie temperature is reduced: 



^ n,,R^n In^ 



\Vo\/Tc 



ln(|l^o|/T.; 



(77) 



A rapidly oscillating exchange interaction Ea (|63|l gives: 



ST,[0] 



An 



"^^°l^'ln(|yo|/T. 



\Vo\/Tc 



(78) 



Tc 3(5+1) 
which is independent of the sign of Vq. 

(3) "Dense" system. 

The most interesting situation is the case of a "dense" 
system of magnetic impurities, when rii ^ (miAi)^/^. 



Since the short-range part of the interaction could have a 
range much larger than the lattice spacing, it need not be 
really dense. This requirement can be rewritten as rii 3> 
(Ai/D)^/^. In a dense system, the wave functions of 
the neighboring impurities overlap strongly, and the main 
exchange contribution arises as a result of this overlap, 
not the RKKY interaction through free carriers. Then 
ferromagnetism arises even when no carriers are present, 
if the Bloembergen-Rowland exchange process dominates 
the physics at short distances {Vq > 0): 



T, 



2S{S+l)n, 



J. 



q=Q 



87r5(5 + l)n, 



VqRI (79) 



Note that Tc does not depend on the carrier concentra- 
tion, and should be much higher than that resulting from 
the RKKY interaction. 

The Bloembergen-Rowland exchange integral may be 
weaker than superexchange. In that case, ferromag- 
netism in a "dense" system is suppressed, and spin glass 
order is favored. This is always the case for indirect band 
gap semiconductors, where the short-range ferromagnetic 
exchange is absent. 

The Bloembergen-Rowland mechanism in case of di- 
rect band gap necessarily leads to an increase of maxi- 
mum Tc as a function of concentration of magnetic im- 
purities rii. 



VII. APPLICATION TO GAAS:MN. 

Application of the Anderson model to a real system, 
such as Gai_a;Mna;As is somewhat more involved than 
the model that we considered above, since Mn ion is in 
3d^ configuration with a spin S = 5/2. This configura- 
tion has 5 d-orbitals. For symmetry reasons, there may 
be more than just one conduction or valence band, which 
is the case for GaMnAs. Then, as we discussed above 
in Section V, one has to take into account all pairwise 
contributions to the exchange integral, Eos. (|28|l . (12911 . In 
particular, there may be unusual contributions to RKKY, 
such as those considered in Section V. In general, one 
has also to sum over all orbitals in the Anderson Hamil- 
tonian, not just spins, and take the Hund's rule, spin 
orbit, and crystal field splitting into account. Different 
orbitals may have different V — s with conduction and 
valence bands. The result of this treatment, however, 
produces the same exchange integrals J{R)', for exam- 
ple, in case when spin-orbit and crystal field splitting is 
neglected, the result will still be given by the exchange in- 
tegral Eqs.l2Hl), 123, with 1^1 1 2 is replaced by l^imP, 
and similarly for IV2P ,where Vim is Vi for m-th orbital. 
The relation between corrections to energy levels and ex- 
change integrals is different in this more realistic case. 
However, the integrals involved are the same, and, as we 
have seen above in section HI, corrections to energy lev- 
els are rather small. Taking Hund's rule into account 
results in replacing S" = 1/2 operator in equations of the 
previous section by S/{2S), with S = 5/2. Thus, this 
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leads to the same results as in the previous section, with 
somewhat redefined V — s. The relationship between the 
energy shift and Jpd will change, but the actual form of 
J(i?) is determined by the energy spectrum only. Instead 
of Vam, we may introduce 

|y„pEE ^-'J^"'' , (80) 

and use the form of exchange integrals that we have ob- 
tained in the previous sections (with S — 5/2). Equiva- 
lently, this would mean that the expression for J{R) in 
terms of Jpd-s (or Jrkky) for conduction and valence 
bands and level position will stay the same as in the pre- 
vious sections. 

For a p-type semiconductor, such as Gai_a;Mna;As, we 
can adopt an inverted picture, where the "empty" band 
is now the valence band (empty of holes) , while the filled 
band is the conduction band. The Anderson Hamilto- 
nian is then easily rewritten in terms of holes. There 
are two types of holes in GaAs - the heavy hole and the 
light hole. Since their masses are very different (O.OSlme 
and 0.5lTOe), the main RKKY exchange contribution is 
produced by the heavy hole. On the other hand, the su- 
perexchange contribution from the light hole band has a 
much larger range, and thus could potentially be more 
important than the superexchange contribution from the 
heavy hole band, or the mixed superexchange contribu- 
tion. However, as we have seen in previous sections, most 
short range contributions for a direct band gap semicon- 
ductor take the form of the second term in Ea. H59|l . al- 
though there may be some variations. The problem is 
that the amplitude of superexchange Vb oc V'^ /D^ cx m? ^ 
where D is the band width! So, for light bands the ef- 
fective range of the interaction is large, but the payback 
is that the amplitude turns out to be small. An easy 
estimate for the light hole band in GaAs shows that for 
V-s of the order of leV this band plays no role in fer- 
romagnetism. The hole mass in the split-off band is 
0.15me. This band could also play an important role 
in the superexchange interaction, although, once again, 
the amplitude for realistic parameters turns out to be 
extremely small. The electron mass is TOge = O.OGStoe 
for the main F-valley The masses and gaps for L- and 



and could change sign as well, at some large doping level. 
If carriers are not present, this exchange integral alone, 
for a > 1, would lead to a saturation or decrease of at 
large doping. On the other hand, when a < 1, ferromag- 
netism gets significantly enhanced at short distances. 
In general, the interplay between various short range 



X- valleys are much larger, so we don't expect them 
to play much role. Thus, we arrive at a simplified pic- 
ture, where only the heavy hole band and F-valley elec- 
trons are relevant. The CFR Mn dP /d^ level, which is 
important for our analysis, is in the conduction band, 
Ai = 1.5eV above the top of the valence band. We can 
see from Eqs (|43f) . (13311 that, for this particular level po- 
sition, the ferromagnetic Bloembergen-Rowland interac- 
tion has the range Rq = Rbr — V\/2mftjjA]", while the 
range of the superexchange is R^^ ~ Rbr/'2.. The am- 
plitude of the Blocmbcrgcn- Rowland term, VsR/Vse — 
V2mge/{2Vimhh), could become comparable or exceed 
the amplitude of heavy hole superexchange. Note that, 
in case of strong short-range ferromagnetic interaction, 
it would be energetically favorable for Mn impurities to 
form ferromagnetic clusters. This, in turn, would re- 
duce the Curie temperature. Clustering of Mn impu- 
rities would make magnetic properties of this material 
crucially dependent on sample preparation. On the other 
hand, antiferromagnetic short-range interactions should 
be stronger at shorter distances, which would potentially 
lead to an exchange integral (in the absence of carriers), 
which changes its sign as a function of the distance be- 
tween impurities. In general, for the particular situation 
when the Mn level is almost at the bottom of conduc- 
tion band, the range of antiferromagnetic superexchange 
is approximately i?o/2, and we may represent the total 
short-range exchange integral in the following form: 

J{R) ^Vo^ h« exp (-2i?/i?o) + exp {^R/Ro)]. (81) 

Here Vq = Vbr > 0, while a ~ m^hVi / {mgeV2) is the 
ratio of superexchange and Bloembergen-Rowland am- 
plitudes (up to a factor of 2), which depends on the 
hybridization of impurity d-level with the valence band 
(V2), of impurity d-level with the heavy hole band, Vhh, 
and the corresponding effective masses. When a > 1 
(which is likely the case here, since ruhh ^ ^ge), the 
exchange becomes antiferromagnetic at short distances 
for R < Rolna. The virial correction to Tc for such ex- 
change integral (assuming Tc is determined mostly by the 
p — d interaction of magnetic impurities and heavy holes) 
is then given by: 



(82) 

I 

contributions leads to a rather complicated physics at 
short distances. While a detailed calculation requires 
precise knowledge of all hybridization parameters from 
the quantum chemistry, we can estimate the Mn concen- 
tration at which the short distance physics becomes im- 
portant by requiring 6Tc/Tc 0.5 in Ea. HVII|l We take 
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FIG. 8: Impurity corrections to Curie temperature. 



the estimate of Jpd ~ IbOeVA^ and Tc ~ llOi^T from 
Ref.[l2, Ai ~ l.beV, and assume that a ~ 1. Then 
Vo = JpamlJiSiT^), and we get: 



VIII. 



2Ro hi (Vo/T,) 



EFFECTS OF DISORDER AND 
INTERACTIONS. 



(83) 



In this section we consider rather briefly effects of dis- 
order and interactions. Since the superexchange and 
the Bloembergen-Rowland exchange interaction are gov- 
erned by high-energy virtual processes, they are indepen- 
dent of disorder. The RKKY interaction, however, gets 
modified. This modification was first considered by de 
Gennes[30l|. who argued that, since the RKKY interac- 
tion at large distances is dominated by the 2pp Kohn 
anomaly wavevector, the vertex corrections are not es- 
sential for the averages over disorder. The long-distance 
power law in the RKKY interaction then gets multiplied 
by an exponential factor exp (— 2i?//), where I is the scat- 
tering length. These effects were indeed taken into ac- 
count by Ohno et al.|l2| in their original paper. Abra- 
hams et al. j3l| , however, have shown, that this is not the 
whole story, since disorder introduces instead a distribu- 
tion of J(i?) at large distances. We note here that the 
long-distance behavior of the RKKY interaction is not 
essential for ferromagnetism. The Curie temperature is 
determined by RKKY exchange at short distances, or 
Jrkky{(1 = 0). Of course, the vertex corrections are es- 
sential for the calculation of the RKKY loop diagram at 
q — 0. Summing all ladder diagrams (0 order in 1/ppl), 
shown in Fig. [SJ leads to a diffuson contribution. 



n(q,w„) = -- 



(84) 



which significantly modifies frequency dependence of 
RKKY at q = 0. Here D is the diffusion coefficient. 

However, the static (w = 0), not dynamic, part of the 
diagram in Fig. [Sldetermines Tc, and it is not changed at 
all. Thus, to the leading order in l/ppl disorder does not 
modify the Curie temperature. The interactions, if not 
too strong, can also be included as the standard Fermi- 
liquid corrections to xa in Eq. Q . Weak localization cor- 



rections (the Cooperon diagrams), however, should mod- 
ify Tc- They can also be included in the same way as the 
standard weak localization corrections to spin suscepti- 
bility (see, for example, Ref-js^), and Eq.Q should still 
be valid. 

Finally, we note that strong exchange interaction JP"* 
could bind holes at Mn sites, forming a shallow (or deep) 
complex magnetic impurity. This effect would reduce 
the hole concentration and the number of free Mn spins, 
and thus lead to a reduction Curie temperature. The 
interactions between these complex magnetic impurities 
would be determined by the overlap of the corresponding 
wavefunctions. 



IX. CONCLUSIONS. 

We have investigated the model of magnetic semicon- 
ductors in which magnetic impurities are treated in the 
framework of the Anderson model. We have shown that 
the effective Hamiltonian of this model is more rich than 
the usual p-d model considered in the literature. Ef- 
fectively, in the Anderson model, wave functions of lo- 
calized impurities develop a "tail", which could be long 
range. When the concentration of magnetic impurities 
is large enough, the overlap of wavefunctions on two dif- 
ferent sites leads to a very strong exchange interaction, 
which is important, and could dominate the physics at 
Mn concentrations as low as 5%. There are two contri- 
butions to this exchange interaction - superexchange, in 
which localized electrons are exchanged trough only one 
type of bands (either conduction or valence bands), and 
the Bloumbergen-Rowland term, when the exchange is 
through both conduction and valence bands. We have 
found that, in case of direct band gap, the Bloemberger- 
Rowland exchange is ferromagnetic. This could lead to a 
dramatic enhancement of Curie temperatures in certain 
magnetic semiconductors with a direct band gap, such as 
GaMnAs. One other important consequence of the indi- 
rect exchange is that, if it is ferromagnetic and reasonably 
long range, doped carriers are not necessary to mediate 
ferromagnetism. This leads potentially to a new class 
of high-temperature magnetic semiconductors, with high 
Curie temperatures determined entirely by the interac- 
tion between localized impurities, not Zener mechanism. 
This emphasizes the effort to search for new materials, 
where ferromagnetism is not carrier-driven (for example, 
driven by the Bloembergen-Rowland mechanism). An- 
other important consequence of the Anderson model is 
that, if there are more than one type of carriers (for ex- 
ample, light and heavy holes in GaAs), the long-range 
RKKY interaction becomes rather complicated, since it 
involves a "mixed" contribution. We have also found 
that, at large doping, the RKKY interaction for the An- 
derson model and the p—d model is different at short dis- 
tances. The effective exchange interaction in the U — oo 
Anderson model for any FS is given by Eas.ltT^ , . 
A numerical solution of this effective Hamiltonian for a 
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given set of parameters (determined from quantum chem- 
istry) would give the answer for Tc{ni) in the most gen- 
eral case. 

We have briefly considered effects of disorder and in- 
teractions. We have shown that when De Gennes[30l| 
approximation ppl ^ 1 is applicable, disorder does not 
modify Curie temperature. Localization effects, however, 
do modify Curie temperature, although their effect could 
be reduced to removal carriers. Finally, in this paper we 
have not considered the effects of mixing of conduction 
and valence bands (such as kp) . These effects should also 
be included in the full description. We should mention 
that application of the Anderson model to GaAs:Mn was 
also considered in Ref. although the limits of the 



Anderson model and their conclusions are different from 
ours. 
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